Polarization-balanced design of heterostructures: Application to AlN/GaN 

double-barrier structures 
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Inversion- and depletion-regions generally form at the interfaces between doped leads (cladding 
layers) and the active region of polar heterostructures like AlN/GaN and other nitride compounds. 
The band bending in the depletion region sets up a barrier which may seriously impede perpendicular 
electronic transport. This may ruin the performance of devices such as quantum-cascade lasers and 
resonant-tunneling diodes. Here we introduce the concepts of polarization balance and polarization- 
balanced designs: A structure is polarization balanced when the applied bias match the voltage 
drop arising from spontaneous and piezeolectric fields. Devices designed to operate at this bias have 
polarization-balanced designs. These concepts offer a systematic approach to avoid the formation 
of depletion regions. As a test case, we consider the design of AlN/GaN double barrier structures 
with AljGai-jN leads. To guide our efforts, we derive a simple relation between the intrinsic 
voltage drop arising from polar effects, average alloy composition of the active region, and the alloy 
concentration of the leads. Polarization-balanced designs secure good filling of the ground state for 
unbiased structures, while for biased structures with efScient emptying of the active structure it 
removes the depletion barriers. 

PACS numbers: 73.21.-b, 73.63.-b, 78.67.-n, 61.46.-w 



I. INTRODUCTION 

Heterostructures of aluminium nitride (AIN) and gal- 
lium nitride (GaN) are key material systems for novel op- 
toelectronic devices.—'^ The large band offset ^ 2eV and 
direct-band gap permits intersubband transitions in the 
near-infrared;^^— which makes them attractive for build- 
ing devices working at fiberoptic wavelengths. They hold 
the promise to extend the operational regime of quan- 
tum cascade (QG) lasers both in the near and the far- 
infraredi^"— The latter due to the large LO-phonon en- 
ergy ujlo ~ 90meV. For wurtzite AlN/GaN structures 
grown in the c-plane direction, the band profile is dis- 
tinguished by huge internal electrical fields. These fields 
arise from the interface charges originating in sponta- 
neous and piezoelectric polarizationii2iH 

AlN/GaN intersubband devices offer great opportu- 
nities, but challenges remain both for growth and de- 
sign. Interface roughness,^'* high dislocation density, and 
scarcity of high-quality substrates limit device perfor- 
mance. Moreover, precise control of multiple quantum 
levels and efficient transport through the structure de- 
mand ultrathin AIN layers. This demand arise because 
the large band offset and heavy conduction electron mass 
~ (0.2 — 0.3)mo weakens the coupling between states 
in adjacent quantum wells. The wide depletion regions 
and two-dimensional electron gases (2DEGs) forming at 
the two sides of an active region , ^1*^1-*^ ^" -^ exacerbates this 
challenge for devices based on perpendicular quantum 
transport. For instance, in resonant-tunneling diodes and 
QG lasers, the depletion regions make it difficult to simul- 
taneously obtain efficient current injection and control 
over the alignment of levels and transitions in the struc- 
ture. Presently, current-driven devices are mostly lim- 
ited to resonant-tunneling diodes, often exhibiting large 
hysteresisJ^ii^"— A related issue exist for structures op- 



erating at zero bias (like absorption structures), espe- 
cially those of few layers: as electrons accumulate in the 
2DEG, the active region is depleted j^'^'^^'i^ 

To tackle some of the design problems, we here in- 
troduce the two related concepts of polarization balance 
and polarization-balanced design. A structure is polar- 
ization balanced when the unscreened field vanish in the 
cladding layers (buffer and cap) layers. We will denote 
these regions as leads if the structure is intended for per- 
pendicular transport. By unscreened we mean the field 
that arise from the polarization-induced interface charges 
(Tpoi and from the applied bias Vbias including dielectric 
screening by valence electrons, but excluding screening 
by conduction electrons and donors. The concept is sim- 
ilar in spirit to strain balance;^ The electrical fields in 
the leads (cladding layers) generally bring about deple- 
tion and inversion regions outside the active region. Since 
the band tailing in the depletion region impedes electron 
transport, polarization balance helps enable efficient in- 
jection into the central active structure. Nonpolar struc- 
tures are always polarization balanced at zero bias. Many 
polar structures are polarization balanced at a certain 
bias. 

Polarization-balanced designs are structures tailored to 
operate at a bias V^^^ resulting in polarization balance. 
In this study, we rely on alloyed leads to retain a large 
design freedom for polarization-balanced designs. A few 
observations substantiated in this paper clarify the con- 
cept and highlight opportunities of such designs. First, 
these designs can, at a desired operational bias V^\^^ ^ 0, 
provide flat-bands in the injection region and in the leads 
in general. There is no depletion barrier. Second, polar- 
ization balance at V^^^^ = is also possible by using a 
suitable choice of composition in the cladding layers. In 
this case, charge accumulates within the active region. 
Third, if charge accumulates in the active region — which 
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Figure 1: Illustration of transport-inhibiting depletion barri- 
ers forming in typical heterostructures comprised of AIN and 
GaN layers. The figure shows conduction-band profile and 
bound states (upper panel), and charge density (lower panel). 
The active regions consist of 1.5 nm AIN barriers sandwiching 
a 2.5 nm GaN layer. Cladding layers consists of doped GaN 
(5 • 10^*cm~'^). The gray background indicates the structure. 
In the upper panel the full curves show the conduction-band 
profile and the dotted indicates the Fermi level. The wavy 
curves display the computed wave function squared of the 
three (quasi-)bound states of the structure. Their offsets are 
given their energy levels. The dashed lines indicate a possi- 
ble profile for an undoped heterostructure. In the lower panel, 
the full curve shows the carrier density, while the dashed lines 
indicate the background doping density. 

can also happen for a finite bias, for instance (in a quasi- 
static picture), when the lowest level of the active region 
falls close to or below the Fermi level — polarization bal- 
ance does not result in flat bands. Fourth, a polar device 
can be operated for a bias before, at, or beyond polar- 
ization balance, that is t/Sbias < "y^b , "y^b ~ </5b , or (^b > V'b 
with (pb = -Vbias and = -^was-'^^ Fifth, nonpolar 
structures cannot operate before, and typically has to 
operate much beyond the bias resulting in polarization 
balance. Sixth, since in polar structures, current-driven 
devices such as quantum-cascade lasers can operate at 
the bias giving polarization balance, with ensuing flat 
bands in the injection region, good injection can be en- 
sured without the need of tailoring graded digital alloys 
in the transport regionJS 

Our calculations of conduction-band profile and 
quantum states presented here are based on the 
Schrodinger-Poisson (SP) equatio n^^i^^ in the effective- 
mass, envelope-function approximation, including non- 
parabolicity and exchange-correlation effects as described 
in Ref . [H 

The plan of the paper is as follows. The second sec- 
tion motivates polarization-balanced designs. We illus- 



trate the formation of depletion and inversion regions for 
an AlN/GaN heterostructure; we demonstrate that the 
depletion barrier seriously inhibit perpendicular current; 
and finally, we take a trial-and-error approach to polar- 
ization balance for AlN/GaN double-barrier structures 
with AlGaN leads of varying alloy concentrations. In the 
third section, we develop theory to guide polarization- 
balanced designs. For an active region consisting of AIN 
and GaN layers strained to Al^Gai^^N leads, a simple 
relation between intrinsic voltage drop, alloy concentra- 
tion i, and average Al composition x in the active re- 
gion is developed. The fourth section contains some sim- 
ple applications: We first demonstrate that choosing a 
polarization-balance design for an unbiased structure se- 
cures good electron filling in the active regions; we then 
turn to a large parameter-space study of a two-level quan- 
tum cascade structure. We find that the use of AlGaN 
instead of GaN leads, significantly enlarge the freedom 
for realizing polarization-balanced designs. The final sec- 
tions hold our discussions, prospects, and conclusions. 

II. DEPLETION BARRIERS IN ALN/GAN 
HETEROSTRUCTURES 

That depletion regions form outside the active re- 
gions in polar heterostructures is well known 
We here observe that these have a huge impact upon 
the function and transport properties of AlN/GaN het- 
erostructures. By calculating the conduction-band pro- 
file of a double barrier we first illustrate the formation of 
such a depletion barrier. Then, by calculating the per- 
pendicular current through a single AIN barrier having 
GaN leads, which allows for a simple and transparent 
analysis, we make evident how severely this barrier in- 
hibits current injection. Finally, we bring to view how 
using AlGaN leads rather than GaN shortens or even re- 
moves these depletion regions. Such an approach opens 
the possibility of avoiding the depletion regions without 
needing huge biases that may prohibit realistic device de- 
sign. This result motivates our theory developments and 
SP calculations in later sections. 



A. Depletion and inversion regions 

Interface charges arise at the boundary between differ- 
ent layers of polar heterostructures as follows: 

a^^i = z ■ (-P'-i + P') = -P'-i + . (1) 

The conduction band profile produced by these charges 
exhibit a characteristic sawtooth-like shape for a layered 
heterostructurei^ i-^^'^'^i-^^'^^i^^i'^" If the lead composition 
and lattice constant is the same on both sides of the 
active region, so that Pq = -P/V) the polar contributions 
to the electrical field vanish in the leads because the net 
charge produced by the polarization-induced charge in 
the active structure is zero: '''poi — 0. 
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The intrinsic voltage drop over the active structure, 
arising from the polarization- induced charges Vpoi, does 
not, in general, vanish. For a structure in equilibrium 
Vhins — some mechanism must compensate for this po- 
tential drop, so that the Fermi level can remain constant 
throughout the structure. For a doped structure, the 
electrical fields rearrange conduction electrons into an 
inversion (accumulation) region on one side of the struc- 
ture, leaving behind a depletion (space-charge) region on 
the other. This charge transfer results in a potential drop 
that cancels the intrinsic voltage drop for a structure 
which is unbiased in equilibrium. If the screening by 
conduction electrons and ionized donors is so weak that 
the Fermi-level falls close to the valence band — as for 
low donor concentration and large potential variations 
— valence electron may tunnel into the conduction band 
creating holes. Hence such polarization-induced doping 
can also compensate for the intrinsic potential drop."^^ 
In our study, the formation of depletion and inversion 
regions in the conduction band is the sole mechanism 
responsible for securing equilibrium. 

Figure [1] illustrates the formation of depletion and 
inversion regions for an unbiased double barrier struc- 
ture, with AIN barriers, and GaN well and cladding lay- 
ers. The conduction-band profile (upper panel) and elec- 
tron density (lower panel) are obtained in a Schrodinger- 
Poisson calculation.^^ We use hard wall boundary con- 
ditions far outside the active region (boundaries not 
displayed) so that a proper bulk develops within the 
cladding layers. The lower panel shows that as electrons 
accumulate on the left, a wide depletion region forms 
(~ 25 nm) on the right. The upper panel shows that 
the resulting band profile exhibits band bending in the 
inversion and depletion regions, and the characteristic 
sawtooth shape in the active part. The ground state 
of the quantum well lies far above (~ leV) the Fermi 
level (dotted lines). In contrast, the energy level of the 
2DEG falls below the Fermi level. In this inversion re- 
gion, several quantised states exist, but electrons accu- 
mulate mostly in the 2DEG of lowest energy .^^i^^ iH The 
dashed lines indicate the conduction-band profile for an 
undoped structure for some arbitrary boundary condi- 
tion. The central role played by the inversion and de- 
pletion regions for this simple heterostructure underlines 
the importance of considering the long-range screening 
as part of a device design. 



B. Tunneling barrier 

The wide depletion barriers seriously reduce tunneling 
current through a heterostructure. A qualitative WKB 
(Wentzel-Kramers-Brillouin) analysis serves as an illus- 
tration. The WKB tunneling probability is as follows: 



Here identifies the classically forbidden regions and the 
angle brackets, an average over the forbidden region. U is 
the potential energy, E, the energy, and m, the effective 
electron mass. The tunneling current is given by the 
Tsu-Esaki formula)^ 
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If Vhias is large, a typical operational condition for 
AlN/GaN structures, electrons are not Pauli blocked 
(which enters through the support function (E) and 
Fermi level /^f)- The tunneling current is then roughly 
proportional to the tunneling probability. We make two 
observations: The heavy mass and large band offset of 
AlN/GaN reduce the tunneling probability compared to 
many other materials. And, a long classically forbid- 
den region reduces the tunneling probability more than 
a high barrier does. In Eq. [5]the term in the exponential 
increases linearly with the length, while the square root 
dampens the dependency on the barrier height. Thus, 
the barrier in the depletion region, which is typically of 
a wide spatial dimension, but of low energy, has a signif- 
icant current-impeding effect. 

The current- voltage curve for a single AIN barrier with 
GaN leads is here calculated using a transmission co- 
effient obtained with a transfer-matrix method for arbi- 
trary barrier shapes;^ and the current / with the Tsu- 
Esaki equation of Eq.[21 This calculation further corrobo- 
rates the WKB argument and attest the extreme current 
sensitivity to the depletion layer width. The AIN bar- 
rier is 1.0 nm wide (two monolayers). A positive back- 
ground charge density of Nd = 5 • lO^^cm"'^ accounts 
for doping of the GaN leads. The conduction-band pro- 
file U{x) is obtained in an effective-mass Schrodinger- 
Poisson description.^^ A stepped Fermi- level induces the 
non-equilibrium charge-density distribution. 

Figure [H displays the current-voltage curves for the 
barrier structure, calculated for conduction-band profiles 
obtained with and without the polar effects. The cur- 
rent increases enormously with bias until about —lAV, 
for which the current in the polar case exceeds that of 
the nonpolar case (which does not exhibit this extreme 
sensitivity). The extreme bias sensitivity has also been 
documented in Refs 15 and 27. The two lower insets 
show the conduction-band profiles in the polar case for 
two sample voltages. It shows that for a voltage of —0.5V 
there is a wide depletion region, while for —1.5 voltage, 
the band is almost flat, with traces of an inversion region 
influencing the right side. It is evident that the reduction 
in sensitivity coincides with achieving flat bands in the 
leads. 

The voltage drop arising from the polarization fields 
is —1.38V. Thus, for an applied bias of Vbias = —1.38V 
the unscreened fields vanish, as the applied bias match 
the voltage drop over the barrier region. For the single 
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Figure 2: Calculated current-voltage curve for a Inm wide 
AIN barrier with GaN leads. The lower green (upper) blue 
curve shows the current calculated with spontaneous and 
piezoelectric polarization included (neglected). The lower in- 
sets (upper inset) show the conduction band profile includ- 
ing (neglecting) these effects at a bias of Vbias = —0.5 and 
— 1.5 V( — 1.5V). The dashed lines identifies the stepped ef- 
fective Fermi-level. 



barrier, the conduction band in the lead region, as ob- 
tained in the SP calculation, is approximately flat when 
the applied bias match the intrinsic voltage drop. 

Where the single barrier in a polar heterostructure 
represents an elementary case of polarization balance, a 
nonpolar heterostructure is the trivial case; Polarization- 
balance is achieved at zero bias. The upper inset in Fig. 
[2] shows that, in the nonpolar case, the depletion and in- 
version regions that arise at nonzero bias do not inhibit 
current as the direction of the slope coincides with the 
movement of the electrons .^^^^ 



C. Biased double barrier structures 

Having studied a single barrier, we next study the 
conduction-band profile of double barrier structures 
which is arguable the simplest structures that exhibits an 
internal quantum structure within the active region. Its 
intended functionality is electroluminescence; the elec- 
tron tunnel resonantly into the upper level of the quan- 
tum well, possibly make a transition to the lower level, 
and exit by tunnelling into the continuum. This toy 
structure allows us to exemplify design features of com- 
plex structures. The few parameters involved in its spec- 
ification also permit brute-force trend studies, such as 
that presented in section 4. 

In designing many different structures, we found that 
a typical feature for AlN/GaN structures with GaN leads 
was that the bias required to overcome the depletion bar- 
riers, and thus get good transport conditions, was so large 
that it resulted in poor confinement of the upper level. 



and very hot electrons in the ejection lead. The left pan- 
els of Fig. [3] illustrates these issues for a double-barrier 
structure. It displays the result of a Schrodinger-Poisson 
calculation for different applied biases, under the assump- 
tion of efficient emptying of the active region; the active 
region and the left lead have the same Fermi level. In 
the cases displayed in panels 'al)' and 'a2)', a depletion 
barrier hinders electron injection. In panel 'a3)', the bias 
is large enough to give flat bands in the leads, and also 
within the quantum well. In this case, states in the in- 
jection lead resonate with the upper level of the active 
region. These energy levels match because the AIN layer 
is appropriately thin for a GaN well of this width. In po- 
lar heterostuctures, the energy levels depend strongly on 
the barrier widths, because of the typically large poten- 
tial drops across the barrier. The 'a3)' structure exhibits 
troublesome features; the potential of the left barrier is 
so low that electrons easily escape prior to making a tran- 
sition, and since the electron is confined in the triangular 
part of the barrier, increasing the well width does not 
help. That the electrons leaving the structure are hot 
can also be a problem. For periodic structures such as 
QC lasers, electrons may never relax before reaching the 
next period. It is also inefficient from a power-usage per- 
spective. 

Using AlGaN leads with a finite Al concentration i, 
instead of GaN leads, reduces the potential drop over 
the active structure at high-transport resonant-tunneling 
conditions since the unscrened field arise from the differ- 
ence in polarization between the layers. As the quan- 
tum transition happens mostly within the active struc- 
ture, the alloy scattering introduced should not severely 
broaden this peak. The middle and right panels, 'bx)' 
and 'cx)', of Fig. |3]show band profiles calculated for al- 
loyed leads. The 'bl)' panel shows that for zero bias and 
X = 0.3, the depletion region is much reduced compared 
to a: = 0.0 in 'al)'. By increasing the bias, as shown 
'b2' panel, the band in the leads eventually becomes flat, 
requiring only about half the bias as with GaN leads as 
shown in 'a3)'. The 'b2)' structure confines the upper 
level much better than if the leads consist of GaN, and 
the electrons tunneling out of the structure remain fairly 
cold. Thus, it outperforms 'a3)'. Both 'b2)' and 'a3)' are 
polarization-balanced designs; they operate close to the 
flat-band condition. 

The 'cl)' structure highlights another important effect: 
Here, the energy of the lower level of the active region 
falls below the conduction band edge of the leads. This 
results in large band filling, with ensuing depletion re- 
gions on both sides of the barrier. Note that for the 'c2)' 
and 'c3)' structures, our assumptions about a stepped 
the Fermi level might be inappropriate. We include the 
structure for the purpose of illustration and complete- 
ness. We find that for this particular active structure, 
the use of Alo.5Gao.5N leads [right panels / 'cx)'] pro- 
hibits flat band in the leads irregardless of applied bias. 

The conduction-band profiles of Fig. |3] demonstrate 
that we can remove the depletion barrier at, and above. 
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Figure 3: Conduction-band profiles and bound states in quantum well of double barrier structures with 0.75 nm wide AIN 



barriers sandwiching an 1.7 nm wide GaN layer for varying bias (fih 
AlGaN leads. 



-Vbiaa and Al concentrations x in the doped (5-10 cm ) 



a certain bias, and that this bias is lowered by using al- 
loyed leads, unless electron accumulate within the active 
structure. Since flat band benefits electron injection, de- 
signing structures to operate at this bias is an excellent 
starting point for making current-driven devices where 
high power is essential. A strategy is therefore to re- 
strict the design space to those structures which are po- 
larization balanced. The next section develops theory for 
conditions that ensures an automatic restriction to such 
polarization-balanced designs. 



III. POLARIZATION BALANCE 

A. Intrinsic voltage drop 

The voltage drop over a heterostructure arising from 
polar effects is found by solving Poisson's equation. The 
solution is 



N 



(5) 



Here Pioad is the polarization of the leads, assumed the 
same in both leads, while Pi and is the polarization 
and dielectric constant of layer i. The polarization of 
a layer has both spontaneous and piezoelectric compo- 
nents: = /^^p ^ ^piczco 

For systems with slowly varying composition or lattice 
constant, an integral formulation is more appropriate: 



/•right lead 

Vpol = dz [P{z) - Plead] — . (6) 

^Icft lead 



Naturally, this expression reduce to Eq. [S] for discrete 
variations of composition and lattice constant. Equa- 
tions [5] or |6] determine the intrinsic voltage drop under 
broad conditions. Polarization balance is obtained when 
it equals the applied bias Vpoi = V^^^. This condition 
generally requires that the two lead compositions are the 
same; otherwise, interface charges does not cancel out, 
and the unscreened fields cannot vanish in the lead re- 
gion. 

If the active region is fully or partially relaxed, the 
lead regions on both sides should be defined so that they 
have the same lattice constant, and thereby piezoelectric 
contributions. 
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B. Heterostructures of two-component alloys: 
AlN/AlGaN/GaN 

The expression for intrinsic voltage drop is here fur- 
ther developed for heterostructures with layers having 
alloy concentrations within the same class of binary al- 
loys. We use AlN/AlGaN/GaN layers as a specific case, 
but the expressions can be applied to other material sys- 
tems. Here, Vegard's law determines the spontaneous 
polarization, P^p{w) = P^g^^i^ — w) + P^^^w, the bulk 
lattice constant, a^uikiw) = (1 — w)aGaN + wqain, as 
well as other basic parameters of Al^^Gai-^jN. The po- 
larization balance concept is not restricted to this ap- 
proximation, but computations are simpler and it allows 
for handy equations to be developed. Vegard's law is, for 
these quantities, a fair approximation for AlGaN layers, 
but it is less appropriate for alloys with larger bowing 
parameters, such as InGaN and InAlN. — 

The piezoelectric polarization can be expressed as 



(7) 



where Xi is the alloy concentration of the layer and yi is 
identified by the bulk alloy concentration corresponding 
to the lattice constant of the layer, that is 



abuik(yi) = at 



(8) 



The heterostructure may have a common lattice constant 
Xi = x' . The piezoelectric material parameters F"^ and 
are 

2D{w) 

[O.A\N — OGaNj 

[a{w) - acaN] 



a(w) 
2D{w) 
a{w) 



(9) 



Here D{w) = 631 (w) - e33{w)ci3{w) / €33(10) with C13 and 
C33 being elastic coefficients and 631 and 633, piezoelectric 
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Figure 4: Piezoelectric material parameters Fa and Fs as a 
function of the alloy concentration of an Al^Gai-^jN layer, 
as given by Eq. O 



coefficients j2S In particular, we identify F^(0) = FcaN, 
F^(0) = 0, andF^(l) = F^(l) = Fain- Figureldisplays 
Ta and Fs as a function of alloy concentration w. While 
F^ is almost linear, F^ exhibits clear bowing. A linear 
approximation for this term is inaccurate. 

With the notation introduced here, we obtain 



N-1 

E 

1=1 



[Pi 



Ff -Plead - 



Pi 



lead 



-rP^'P 
■^-^AIN 



(i-s)^GaN + yr^(^)-r^(5) 



(10) 



Here x is the alloy concentration of the leads and y is the 
bulk AlGaN concentration corresponding to the lattice 
constant of the leads (Eq. |8]) . 

Two consistency checks of Eq. [TU]are readily available: 
First, if a layer (for instance the leads) is relaxed, Xi = yi, 
its piezo-electric component, Ppiczoixi, Xi) — Xir^{xi) — 
T^{xi) = 0, vanish. Second, the intrinsic voltage drop 
does not depend on the starting point of the lead region: 
If we define part of a lead (layer k) as a layer within the 
active region, the intrinsic potential drop V^oi remains 
the same, because the additional k term yields zero: Pk — 

Plead = 0. 



C. AlN/GaN layers strained to AlGaN leads 

If the active region is comprised of AIN and GaN layers, 
and these are strained to the leads (x' = x = y), there is 
a simple relation between the lead alloy concentration i, 
intrinsic voltage drop over active region length Vpoi/L, 
and average concentration in the active region: 



x = Y^LfVL. 

The relation can be stated as follows: 
~ _ ~Vpoi/L + Ax 



Ax + B{1 - x) 



The material specific values A and B are given by 



A^{P^f^~P, 



jGaN 
sp 



1 



B 



P 



AlN 



p: 



GaN 



F 



GaN' 



EAIN 
I 

CGaN 



(11) 

(12) 

(13) 
(14) 



with materials-specific values listed in Table HI In deriv- 
ing the relation (|12p between intrinsic voltage drop and 
lead concentration, we have used the Vergard's law as in 
the previous subsection. 

Figure [5] displays the relation (|12p between the average 
composition x and lead allow concentration x for different 
Vpo\/L. By setting the applied bias according to these 
contours, we get polarization balance, that is, we remove 
the unscreened fields in the leads. 

We confirm the consistency of Eq. [T^ by showing that 
two polarization balanced active regions are also natu- 
rally described as single active region. We consider two 
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Table I: Polarization-related material quantities introduced 
here, based on the parameters of Ref. [39l . 



A [V/A] 


B [V/A] 


Tain [C/m^] 


TcaN [C/m^] 


-0.138 


-0.094 


0.048 


0.025 




LO 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 
Average active structure composition x 



Figure 5: The intrinsic voltage drop per length Vpoi/ L [V/nm] 
of the active region as a function of the average composition 
X of the active region and the alloy concentration of leads x. 
The active region is strained to the leads. 



active regions with average compositions xi and X2, of 
length Li and L2 producing respectably a intrinsic volt- 
age drop of Vi , V2 for the same lead alloy concentration 
X. We test whether x is the solution to the following 
equation (with unknown z) describing the two regions in 
terms of a single active region: 



V1+V2 



A 



XiLi + X2L2 



L1+L2 
■^a;iii + X2L2 



L1+L2 
+ B{1- 



xiLi + X2L2 



L1+L2 Li+L2 

Singling out the Li and L2 terms, we obtain 

Li [-V1/L1 + Axi-z {Axi + B{1 - xi))] - 
L2 [-V2/L2 + Ax2~z {Ax2 + B{1 - X2))] = 



(15) 



0. (16) 



Since z = 5; is a solution, the two descriptions are equiv- 
alent. 




Figure 6: Conduction-band profile and bound states of an 
optically active quantum well with a polarization-balanced 
design for zero bias. It consists of 2.5 nm wide AIN barriers, 
a 1.5 nm wide GaN well layer and cladding layers consisting of 
Alo.63Gao.37N. Cladding layers are doped to 5-10^*cm~^. The 
light gray [green] (white, gray) background indicate AlGaN 
(GaN,AlN). The dashed [purple] lines show the band profile 
of an undoped structure. The wavy curves indicate the wave 
function squared of the bound states offseted by their energy. 
The dotted line represents the Fermi level. 



GaN quantum well in-between AIN barriers and AlGaN 
cladding layers with alloy concentraction chosen for po- 
larization balance at zero bias. There is a significant band 
filling of the active region which generates depletion re- 
gions on both sides of the active region. This contrasts 
the depletion and inversion regions outside the depleted 
active region for GaN cladding layers, as displayed in 
Fig. [T] and the 'la)' panel of Fig. [31 The dashed lines 
confirms the polarization balance equations, giving fiat 
band in the cladding layers for an undoped structure. 

The elimination of unscreened fields in the cladding 
layers (polarization balance) cause charging of the active 
region for unbiased structures because the energy of the 
lowest quantum well state falls below the condition band 
edge. The band bending in the cladding layers under- 
line that merely fulfilling polarization balance does not 
necessarily result in a fiat band. In the unbiased case po- 
larization balance instead ensures efficient filling of the 
active region, rather than of the 2DEG at the interface. 
This result relates to an earlier observation that using 
the same alloy composition in the active region as in the 
cladding layers achieves good control over the electron 
population<^iil 



IV. EXAMPLES OF APPLICATIONS 

A. Double barrier structure polarization balanced 
for zero bias 

Figure H] displays the conduction-band profile and 
bound states of an unbiased structure comprised of a 



B. Design of quantum-cascade structure 

If an active region is approximately charge neutral 
at operating conditions, the polarization-balanced design 
succeeds to remove depletion and inversion regions. This 
could be the case for heterostructures with efficient emp- 
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Figure 7: (color) Design criteria, transition energy, and applied bias contours for polarization balanced two-level QC structure 
as a function of well and barrier widths for different AlsGai-jN alloy concentrations of the leads: x — 0, 0.2, 0.3, and 0.4. 
The numbers along the dashed curves mark the absolute value of applied bias (Vbias [V]), set to achieve polarization-balance 
(using Eg. I12p . while those on the full, mark the transition energy AE12 [meV] between the two levels in the well. The white 
areas indicate where all design criteria are met; the shaded semi-transparent where they are not: The blue areas indicate where 
confinement of the upper level is poor. The yellow, where ejection of the lower level is problematic. In the green (light green) 
areas, the injection is problematic because the energy of the upper level falls below (far above) the band edge of the injection 
lead. The inset shows the band profile and the absolute-square of the localized wave-functions offset by their energy for the 
simple QC structure, and the shaded regions indicates design criteria. 



tying of the active region. By adjusting the bias over a 
structure, we can always make a structure polarization 
balanced; however, the required bias can be unsuited for 
meaningful designs of the active structure. This was illus- 
trated for GaN leads in our trial-and-error study detailed 
in Fig. 13 

The additional freedom of design gained by using al- 
loyed leads for polarization balanced structures is here 
demonstrated with a large parameter study of a simple 
two-level quantum-cascade (QC) structure. The active 
region consists of two AIN barriers sandwiching a GaN 
well as in previously studied structures. We vary well 
(Lw) and barrier (Lbar) widths and lead alloy concentra- 
tion X. The bias is always adjusted to obtain polarization 
balance. In calculating energies and wave functions, we 
only consider the active region. This is since we found 
only a tiny difference between the band profile as ob- 
tained with and without doping of the leads, given that 
charging of the active region is avoided. 

The merits of a design depend on the intended ap- 
plication; for instance, gain and threshold current are 
key quantities for lasers. Rather than considering an 
application-specific assessment of our QC structure, we 
consider three intuitive, but somewhat crude, design cri- 



teria. The first is injection of the upper level: the energy 
should lie above the conduction band edge of the injec- 
tion lead, but no more than 200 meV above. This upper 
bound is somewhat arbitrary; in practice it depends on 
the band filling of the leads and level broadening. How- 
ever, the applied bias can be adjusted to accommodate 
slight misalignment. The second is confinement of the 
upper level, the energy should fall below the triangular 
portion of the barrier to increase the fraction of the elec- 
trons making a transition to the lower level by avoiding 
the faster Fowler-Nordheim tunneling. The third is ejec- 
tion of the lower level: the energy should lie above the 
conduction band edge of the exit lead. Efficient ejection 
is essential to avoid charge-buildup causing depletion re- 
gions in the injection lead and Pauli blocking of optical 
transitions. Finally, we aim to design a structure oper- 
ating at the fiberoptic frequency ^ SOOmeV. 

Figure [7] details the results of our large parameter 
study of the simple QC structure. It shows that the 
structures with GaN leads does not fulfil the design cri- 
teria for a transition energy of 800 meV. The bias over 
the structure is also generally about twice the transition 
energy causing ultrahot electrons to exist in and around 
the structure. For larger alloy concentration in the leads. 
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the acceptable parameter window moves to wider barri- 
ers and thinner wells. For x = 0.2, and x = 0.3, all 
design targets can be accommodated for the fiberoptic 
frequency. Increasing lead alloy concentrations further, 
the allowed window contracts, and for x = 0.5 no param- 
eters supports all design criteria. 



V. DISCUSSION 
A. Conduction-band regimes 




Figure 8: Overview of considered structures: Different 
conduction-band regimes as function of applied voltage 
(horizontal-axis) and an abstract structure space (vertical- 
axis). The red curve marks polarization balance conditions. 
For structures along along this curve, the leads bands are 
flat unless charge accumulates within the active region, which 
happens where the curve enters the gray (orange) region. 
The flgure number of structures considered in this study are 
marked on various spots. All cases considered in Fig. 7 are po- 
larization balanced. Those within the dashed circles fulfil cer- 
tain chosen design conditions. This makes them polarization- 
balanced designs. 

Figure |8] provides an overview of the different struc- 
tures and bias conditions considered in this study and 
attempts to sort these into different conduction-band 
regimes. The horizontal-axis indicate different applied 
bias over the structure. The vertical-axis represents an 
abstract structure space; for instance, increasing the Al 
concentration of the cladding layer for a given active 
structure corresponds to moving downwards. The (blue) 
curve marks a structure in polarization balance. The 
gray (orange) region indicates that electrons accumulates 
within the active region, with ensuing depletion barriers 
on both sides of the active region. 

We identify the position of structures studied here: In 
Fig. Ql we considered a structure far from polarization 
balance and with an empty active region. It exhibited a 
wide depletion region on the right side of the structure 
and a well-filled inversion region on the left. The panels 



of Fig. displayed three different structures, under three 
different voltage conditions. The structures under vary- 
ing bias is indicated by the three semi-transparent axes. 
Two of the panels ['a3'),'b2)'] showed structures under 
a bias that made them approximately polarization bal- 
anced. In two of the structures much charge accumulates 
within the active region ['cl)','c2)'] as the Fermi-level was 
above the lowest level of the active region. The structure 
showed in Fig. [^was designed for polarization balance at 
zero bias according to the theory developed in section 3. 
This secures excellent charge accumulation in the active 
region. Fig. displayed key results of numerous struc- 
tures at a bias giving polarization balance. In Fig. [3 all 
these results scatter along the curve indicating polariza- 
tion balance. Some of the designs resulted in filling of 
the active region (placing them within the gray region) 
and are therefore among the failed designs. The purpose 
of the dashed circles is to indicate (in an abstract sense) 
successful polarization-balanced designs according to the 
selected design criteria. 

The placement of the results of the single barrier dis- 
played in Fig. [2] are not shown here. We note that since 
the single barrier does not have any internal structure, 
accumulation within the active structure does not oc- 
cur. A single barrier can not be designed for polariza- 
tion balance at zero bias. The nonpolar case is always 
polarization-balanced at zero bias, but for such struc- 
tures, polarization balance does not guarantee charge ac- 
cumulation within the active region. 



B. Prospects and challenges 

The theory developed for polarization balance may 
be useful for designing structures not intended as 
polarization-balanced designs, since it can also be used 
to identify the flat-band condition. Designing a device 
to operate at an applied bias larger than that resulting 
in polarization balance may be advantages for some pur- 
poses. For instance, to ensure an accumulation layers in 
the injection lead, making its operation more like that 
of a nonpolar structure.'*'^ Nevertheless, if an effective 
continuum should arise in the injection region, the bias 
across it should be small, to avoid the formation of a 
Stark-ladder4M^ 

The huge strain in the structures considered in this 
study might render them unsuitable for actual applica- 
tions, although there are numerous of attempts of making 
resonant-tunneling diodes using GaN leads which have 
even larger strains.i^"— i^i^^ Using AlGaInN layers could 
be helpful, but complicates growth. With this additonal 
freedom, one could design structures that are both polar- 
ization and strain balanced.^'* Only in special accidental 
cases can one fulfil both these criteria using layers within 
the same class of binary alloys (such as AlGaN) . 

Lattice matching by using AlInN barriers and GaInN 
wells"^^ removes strain completely. However, with only 
two different alloy compositions of the layers the generic 
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polarization balance approach is impossible, since it in 
general requires layers of at least three different compo- 
sitions. We therefore mention that polarization balance 
concepts can also be developed for leads that are super- 
lattices. Such an approach is in the spirit of digital alloys. 
It would allow for structures that are both polarization 
balanced and lattice-matched (or the more relaxed strain- 
balance condition), for instance by combining Al82lni8N 
and GaN4ii^ Since strain is zero, we only balance the 
spontaneous part of the polarization. 

Aside from presenting a solution to the current in- 
hibiting effect of band tailing in the depletion region, 
polarization-balanced designs present an additional ad- 
vantage: The shape of the conduction band does not de- 
pend on inversion and depletion regions. This has the 
consequence that the conduction-band profile essentially 
decouples from the specific doping and carrier concen- 
tration. This reduces the importance of an exact control 
over the donor density for the functionality of a device, 
which is also helpful since the number of ionized donors 
is complicated to calculate in heterostructures, or when 
the binding energy is similar to the thermal energy4^ 
We note that in our quasi-static study, effects like charge 
trapping and formation of field domainsii are not in- 
cluded, and these could complicate design considerations. 

For emission structures of short optical wavelength, 
and in particular for QC lasers made up of an alternating 
series of optically active and interjecting current injec- 
tions regions, the possibility to directly engineer zero-bias 
conditions in the injection region (or in the leads) points 
to an additional strong advantage of polar heterostruc- 
tures with polarization-balanced designs. In contrast to 
conventional nonpolar structures, it is not necessary to 
grade the injection layers to compensate for the field. A 
large applied bias does not inevitably result in a large 
electrical fields in the leads or, for that sake, in all of the 
active (transport-limiting) region. In a good design, the 
potential drop can be engineered to occur primarily in 



the region where the optical transition takes place. We 
conjecture that this advantage can result in a significant 
simplification of the tailoring of the quantum wells and 
barriers that ensures efficient injection in devices like QC 
lasers. 



VI. CONCLUSION 

The polarization balance concept represents a guiding 
principle for robust designs of wurtzite AlN/GaN het- 
erostructures grown in the c-direction. For high-power 
current-driven devices, it helps eliminate the depletion 
barrier impeding the perpendicular current. It also sim- 
plifies design since one can, to a larger extent, focus on 
the active structure. This is because the infiuence of the 
doping levels on the conduction-band profile is substan- 
tially reduced as the conduction band is flat in the leads. 
Furthermore, the flat bands in the leads makes grading 
of injection regions in QC lasers unnecessary. For such 
designs of unbiased structures it instead secures a large 
electron population in the active region. 

Our study has been restricted to AlGaN, since the large 
mass and band offset accentuates the current inhibiting 
effect of the depletion regions. However, the concept nat- 
urally generalize, not only to nitrides^ including InN and 
BN and their alloys, but also to heterostructures of other 
polar materials such as ZnMgO/Zn04i 
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